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•  Properties  of  Nonholonomic  Systems: 

•  Energy  Conservation 

Hamiltonian:  Yes.  Nonholonomic:  Yes. 

•  Momentum  Conservation 

Hamiltonian:  Yes,  Noether’s  Theorem.  Nonholonomic:  No,  Momentum  Equation 

•  Measure  (volume  )  Preservation 

Hamiltonian:  Yes.  Nonholonomic:  No,  in  general 

•  Stability 

Hamiltonian:  Never  asymptotic.  Nonholonomic:  Can  be  asymptotic. 

•  Key:  Almost  Poisson  structure,  Nonvariational 
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•Nonholonomic  Equations  of  Motion 

See  e.g.  Bloch,  Krishnaprsad,  Marsden  and  Murray  [1996]  and  Zenkov,  Bloch  and  Marsden 
[1998],  Bloch  and  Crouch  [1995]  and  other  references  in  these  papers. 

•The  Lagrange-d’Alembert  Principle 

•  Consider  a  system  with  a  configuration  space  Q.  local  coordinates  q'  and  m  nonintegrable 
constraints 

sa  +  Aaa (r,  s)ra  =  0 

where  q  =  (r,  s)  G  Mn_p  x  which  we  write  as  ql  =  (ra,  s°),  where  1  <  a  <  n  —  p  and 
1  <  a  <  p. 

•  Lagrangian  L(q\  q1). 

Equations  of  motion  given  by  Lagrange-d’Alembert  principle. 

Definition  0.1  The  Lagrange-d’Alembert  equations  of  motion  for  the  system  are 
those  determined  by 

S  I  L(q\  ql)  dt  =  0, 

J  a 

where  we  choose  variations  Sq(t)  of  the  curve  q{t)  that  satisfy  Sq(a)  =  Sq(b)  =  0  and  5q(t) 
satisfies  the  constraints  for  each  t  where  a  <  t  <  b. 


•  This  principle  is  supplemented  by  the  condition  that  the  curve  itself  satisfies  the  constraints. 


•  Note  that  we  take  the  variation  before  imposing  the  constraints;  that  is,  we  do  not  impose 
the  constraints  on  the  family  of  curves  defining  the  variation. 

•  Equivalent  to: 


-SL=  Sg'  =  0 

\dtdql  dql  J 


for  all  variations  Sql  =  (hra,  5sa)  satisfying  the  constraints  at  each  point  of  the  underlying 
curve  q(t ),  i.e.  such  that  Ssa  +  Aaa8ra  =  0. 

Substituting: 


(. 


d  dL  dL 


\dtdra  dr01 

for  all  a  =  1, . . . ,  n  —  p. 

Combined  with  the  constraint  equations 


=  k( 


d  dL  dL 


\dtdsa  dsc 


for  all  a  =  give  the  complete  equations  of  motion  of  the  system. 


Useful  way  of  reformulating  equations  (0.2)  is  to  define  a  constrained  Lagrangian  by  substi¬ 
tuting  the  constraints  (0.3)  into  the  Lagrangian: 

Lc(ra ,  5a,  fa)  :=  L(r“,  5a,  r“,  -A“(r,  s)fa ). 


The  equations  of  motion  can  be  written  in  terms  of  the  constrained  Lagrangian  in  the  following 
way,  as  a  direct  coordinate  calculation  shows: 


ddLc 
dt  dra 


dLc 

dra 


+  Aa  dsa 


dL_ 

dsb 


na(5r 


i 


where  Bba  >  is  dehned  by 


•  The  Falling  Rolling  Disk 


This  is  a  system  which  exhibits  stability  but  not  asymptotic  stability. 

Denote  mass,  the  radius,  and  the  moments  of  inertia  of  the  disk  by  m,  i?,  A,  B. 


L  = 


m 


+ 


(£  —  R((j>  sin  0  +  ip))z  +  gz  sin  6  +  (77  cos  6  +  R6\ 
1 


2  L 


A(6Z  +  <pz  cos2  9)  +  B((j>  sin  9  +  ip)  —  mgR  cos  9, 


where  £  =  x  cos  (p  +  y  sin  (p  +  Rip  and  g  =  —  x  sin  (p  +  y  cos  (p,  while  the  constraints  are  given 
by 


x  =  —ipR  cos  <p1 


y  =  —  ^Rsin0. 


Other  systems: 


Figure  0.2:  The  Chaplygin  sleigh  is  a  rigid  body  moving  on  two  sliding  posts  and  one  knife  edge. 


Figure  0.4:  The  rattleback. 


•The  Chaplygin  Sleigh 

•  Perhaps  the  simplest  mechanical  system  which  illustrates  the  possible  dissipative  nature  of 
energy  preserving  nonholonomic  systems. 

Compare  the  sleigh  equations  to  the  Toda  lattice  equations. 


Figure  0.5:  The  Chaplygin  sleigh  is  a  rigid  body  moving  on  two  sliding  posts  and  one  knife  edge. 

Equations: 


v  =  auj2 

ma 2 

Cj  =  - - 

1  +  maz 

Equations  have  a  family  of  relative  equilibria  given  by  (u, u)\v  =  const,  u  =  0. 

Linearizing  about  any  of  these  equilibria  one  finds  one  zero  eigenvalue  and  one  negative 


eigenvalue. 

In  fact  the  solution  curves  are  ellipses  in  v  —  uj  plane  with  the  positive  u-axis  attracting  all 
solutions. 
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Figure  0.6:  Chaplygin  Sleigh/2d  Toda  phase  portrait. 


Normalizing,  we  have  the  equations 

v  =  (jJ2 
Cj  —  —vuj  . 

Scaling  time  by  a  factor  of  two  have:  Chaplygin  sleigh  equations  are  equivalent  to  the  two- 
dimensional  Toda  lattice  equations  except  for  the  fact  that  there  is  no  sign  restriction  on  the 


variable  uj.  Hence  can  be  written  in  Lax  pair  form  and  solved  by  the  method  of  factorization. 


•The  Toda  Lattice 

Interacting  particles  on  the  line. 

Non-periodic  finite  Toda  lattice  as  analyzed  by  Moser  [1974]: 


Hamiltonian  equations: 


H(x,y) 


n  7i—l 

k— 1  k— 1 


yk 


dH 

o  —  ///•■ 

oyk 

c)H 

_ _ _  exk-i~xk 


g^/c  xk—  1 


where  assume  e^0  Xl  =  eXn  Xn+1  =  0. 
Flaschka: 


Get: 

&k  ^A:(^/c+l  3  k  I3  '  '  '  3  kl 

h  =  2 (af  -  aPi) ,  k  =  1,  •  •  •  ,n 
with  the  boundary  conditions  ao  =  an  =  0  and  where  the  a\  >  0. 


Matrix  form: 


-rL  =  \B,L]  =  BL-LB, 

dt  1  1 

If  N  is  the  matrix  diag[l,  2,  •  •  •  ,  n\  the  Toda  flow  can  be  written 


L  =  [L,  [L,  N]_  . 


Shows  flow  also  gradient  (on  a  level  set  of  its  integrals). 

•  Double  bracket  form  of  Brockett  [1988]  (see  Bloch  [1990],  Bloch  Brockett  and  Ratiu  [1990, 
1992]). 

•The  Two-dimensional  Toda  Lattice 

In  two-dimensional  case  matrices  in  the  Lax  pair  are 


L 


(  0  a\  \ 

V  ~ai  °  ) 


Equations  of  motion: 


b\  =  2  a\ 
d\  =  —2a\bi 

For  initial  data  b\  =  0,  a\  =  c,  explicitly  carrying  out  the  factorization  yields  explicit  solution 


sinh  2  ct 
cosh  2  ct 


a\(t)  = 


cosh  2  ct 


Symmetries  Symmetries  play  an  important  role  in  onr  analysis.  Suppose  we  are  given  a 
nonholonomic  system  with  Lagrangian  L  :  TQ  — >•  M,  and  a  (nonintegrable)  constraint  distri¬ 
bution  V.  We  can  then  look  for  a  group  G  that  acts  freely  and  properly  on  the  configuration 
space  Q.  It  induces  an  action  on  the  tangent  space  TQ  and  so  it  makes  sense  to  ask  that  the 
Lagrangian  L  be  invariant.  Also,  one  can  ask  that  the  distribution  be  invariant  in  the  sense 
that  the  action  by  a  group  element  g  E  G  maps  the  distribution  Vq  at  the  point  q  £  Q  to  the 
distribution  Vgq  at  the  point  gq.  If  these  properties  hold,  we  say  that  G  is  a  symmetry  group. 
The  manifold  Q/G  is  called  the  shape  space  of  the  system  and  the  configuration  space  has  the 
structure  of  a  principal  fiber  bundle  7 t  :  Q  —>  Q/G. 

Geometry  of  Nonholonomic  Systems  with  Symmetry 

The  group  orbit  through  a  point  q1  an  (immersed)  submanifold,  is  denoted 

Orb(g)  :=  {gq  |  g  €  G}. 

Let  g  denote  the  Lie  algebra  of  the  Lie  group  G.  For  an  element  £  E  g,  we  denote  by  £g  the 
vector  field  on  Q  arising  from  the  corresponding  infinitesimal  generator  of  the  group  action, 
so  these  are  also  the  tangent  spaces  to  the  group  orbits.  Define,  for  each  q  E  Q,  the  vector 
subspace  gg  to  be  the  set  of  Lie  algebra  elements  in  g  whose  infinitesimal  generators  evaluated 
at  q  fie  in  both  Vq  and  Tq(Orb(q))\ 

09  :=  {£  e  0  I  63(g)  e  Vq  n  r„(Orb(g))}  . 

The  corresponding  bundle  over  Q  whose  fiber  at  the  point  q  is  given  by  g-i,  is  denoted  by  g  . 


Reduced  dynamics.  Assuming  that  the  Lagrangian  and  the  constraint  distribution  are 
G-invariant,  we  can  form  the  reduced  velocity  phase  space  TQ/G  and  the  reduced  constraint 
space  V / G.  The  Lagrangian  L  induces  well  defined  functions,  the  reduced  Lagrangian 

l  :  TQ/G  ->  M 

and  the  constrained  reduced  Lagrangian 

lc  :  V/G  ->  K, 

satisfying  L  =  l  o  7 ttq  and  L \x>  =  lc  0  where  7 ttq  :  TQ  — >•  TQ/G  and  nz>  :  P  — >■  V/G 
are  the  projections.  By  general  considerations,  the  L agr angc- d ' Alernb crt,  equations  induce  well 
defined  reduced  equations  on  V /G.  That  is,  the  vector  field  on  the  manifold  V  determined  by 
the  Lagrange-d’ Alembert  equations  (including  the  constraints)  is  G-invariant,  and  so  defines 
a  reduced  vector  field  on  the  quotient  manifold  V/G.  Call  these  equations  the  Lagrange- 
d' 'Alembert- Poincare  equations. 


Let  a  local  trivialization  be  chosen  on  the  principle  bundle  i t  \  Q  Q/G,  with  a  local 
representation  having  components  denoted  (r,  g).  Let  r,  an  element  of  shape  space  Q/G ,  have 
coordinates  denoted  ra,  and  let  g  be  group  variables  for  the  fiber,  G.  In  such  a  representation, 
the  action  of  G  is  the  left  action  of  G  on  the  second  factor.  The  coordinates  (r,  g)  induce  the 
coordinates  (r,  r,  £)  on  TQ/G ,  where  £  =  g~lg.  The  Lagrangian  L  is  invariant  under  the  left 
action  of  G  and  so  it  depends  on  g  and  g  only  through  the  combination  £  =  g~lg.  Thus  the 
reduced  Lagrangian  l  is  given  by 

l(r,r,0  =  L(r,g,r,g). 

Therefore  the  full  system  of  equations  of  motion  consists  of  the  following  two  groups: 

1.  The  Lagrange-d’ Alembert- Poincare  equation  on  V /G  (see  theorem  0.2). 

2.  The  reconstruction  equation 


g  =  gC 


The  nonholonomic  momentum  in  body  representation.  Choose  a  g-dependent 
basis  eA{q)  for  the  Lie  algebra  such  that  the  first  m  elements  span  the  subspace  Qq  in  the 
following  way.  First,  one  chooses,  for  each  r,  such  a  basis  at  the  identity  element  g  =  Id,  say 

ei(r),  e2(r), . . . ,  em(r),  em+1  (r), . . . ,  ek(r). 

Now  define  the  body  fixed  basis  by 

eA(r,g)  =  Ad^  eA(r). 

Then  the  first  m  elements  will  indeed  span  the  subspace  Qq  since  the  distribution  is  invariant. 
We  denote  the  structure  constants  of  the  Lie  algebra  relative  to  this  basis  by  CAB. 

To  avoid  confusion,  we  make  the  following  index  conventions: 

1.  The  first  batch  of  indices  range  from  1  to  m  corresponding  to  the  symmetry  directions  along 
constraint  space.  These  indices  will  be  denoted  a,  6,  c, . . . . 

2.  The  second  batch  of  indices  range  from  m+1  to  k  corresponding  to  the  symmetry  directions 
not  aligned  with  the  constraints.  Indices  for  this  range  will  be  denoted  by  a/,  b' ,  <d, . . . . 

3.  The  indices  A,  B,  C, . . .  on  the  Lie  algebra  g  range  from  1  to  k. 

4.  The  indices  a:  /?, . . .  on  the  shape  variables  r  range  from  1  to  a.  Thus,  a  is  the  dimension 
of  the  shape  space  Q/G  and  so  a  —  n  —  k. 

The  summation  convention  for  all  of  these  indices  will  be  understood. 


Assume  that  the  Lagrangian  has  the  form  of  kinetic  minus  potential  energy,  and  that  the 
constraints  and  the  orbit  directions  span  the  entire  tangent  space  to  the  configuration  space: 

Vq  +  r9(Orb(?))  =  TqQ. 

Then  it  is  possible  to  introduce  a  new  Lie  algebra  variable  Q  called  the  body  angular  velocity 
such  that: 

1.  Q  =  *4rT£ ,  where  the  Lie  algebra  valued  form  A  =  AAeA(r)dra  is  called  the  nonholonomic 
connection  (see  Bloch  et  al.  [1996]  for  details). 

2.  The  constraints  are  given  by  Q  E  span{ei(r), . . . ,  em(r)}  or  Qm+1  =  •  •  •  =  Qk  =  0. 

3.  The  reduced  Lagrangian  in  the  variables  (r,  r,  Q)  becomes 

l{ra,  ra ,  ttA)  =  ^ gayrar p  T  ab^A^B  +  K'Aatta'  -  U(r).  (0.8) 

Here  gay  are  coefficients  of  the  kinetic  energy  metric  induced  on  the  manifold  Q/G ,  Iac  are 
components  of  the  locked  inertia  tensor  defined  by 


(i(OC  n)  =  {{fam)), 

where  {(• ,  ■))  is  the  kinetic  energy  metric.  The  coefficients  A ara  are  defined  by 

d2l  d2l 


^ d(,a'dr 


nf. 


a 


dO'di 


B  ° 


The  constrained  reduced  Lagrangian  becomes  especially  simple  in  the  variables  (r,  r,  fi): 

lc  =  \gaprarp  +  \lab^b  ~  u.  (0.9) 

We  remark  that  this  choice  of  Q  block-diagonalizes  the  kinetic  energy  metric,  i.e.,  eliminates 
the  terms  proportional  to  Qara  in  (0.9). 


The  nonholonomic  momentum  in  body  representation  is  defined  by 


dl  _  dlc 
dtta  ~  dQa’ 


a  =  1, . . . ,  m. 


Notice  that  the  nonholonomic  momentum  may  be  viewed  as  a  collection  of  components  of  the 
ordinary  momentum  map  along  the  constraint  directions. 

The  Lagrange-d’Alembert-Poincare  equations.  As  in  Bloch  et  al.  [1996],  the 
reduced  equations  of  motion  are  given  by  the  next  theorem. 


Theorem  0.2  The  following  reduced  nonholonomic  Lagrange-d’ Alembert- Poin¬ 
care  equations  hold  for  each  1  <  a  <  a  and  1  <  b  <  m: 


Here  and  below  lc(ra1  ra ,  Qa)  is  the  constrained  Lagrangian,  and  Ibd  and  7ay  are  the  inverse  of 
the  tensors  I|0<?  and  respectively.  We  stress  that  in  general  Ibd  Ibd  and  7a/c/  I a/c/. 


The  coefficients  uaJ^  T>(ha ,  T>iJap1  JCa^  are  given  by  the  formulae 

f)  AC  P)  AC 

nc  _  UJXa  nC  aAaB.C  aA  C  aA 

Qr(3  Qra  ° BA^a^f3  +  7 Afi^a  fAo^fj  j 


vcba  =  ~(ccAb  -  cibiMr'c)A£  +  c<ab \daiac  +  7L  -  YbaBA 

'DbaP  =  ^c'P^lba  ~  CCaiA. 

! 


a  c 


JAb^t)i 


fca/3^  ~  ^c'^af3-> 


0.12) 


and  the  coefficients  are  defined  by 

dcjj  q 

Qra  —  ffi baeC • 

Equations  (0.10)  and  (0.11)  generalize  the  equations  of  motion  in  the  orthogonal  body  frame 
(see  Bloch  et  al.  [1996]).  Here  we  no  longer  assume  that  the  body  frame  is  orthogonal. 


•Almost  Poisson  Systems 

Recall: 

Definition  0.3  An  almost  Poisson  manifold  is  a  pair  (M,  {, })  where  M  is  a  smooth 
manifold  and  (i){, }  defines  an  almost  Lie  algebra  structure  on  the  C°°  functions  on  M , 
i.e.  the  bracket  satisfies  all  conditions  for  a  Lie  algebra  except  that  the  Jacobi  identity  is 
not  satisfied  and  (ii)  {, }  is  a  derivation  in  each  factor. 

If  in  addition  Jacobi  satisfied,  Poisson  manifold. 

An  almost  Poisson  structure  on  M  will  be  Poisson  if  its  Jacobiator ,  defined  by 

J(f,g,h )  =  {{f,g},h}  +  {{g,h},  f}  +  {{h,  f}g} 

vanishes. 

•  One  can  define  an  almost  Poisson  vector  field  on  M  by 


•  “Hamiltonian”  Formulation  of  Nonholonomic  Systems 

Nonholonomic  systems  are  almost  Poisson. 

Start  on  the  Lagrangian  side  with  a  configuration  space  Q  and  a  Lagrangian  L  (possibly  of 
the  form  kinetic  energy  minus  potential  energy,  i.e., 

1 

L(q,q)  =  W FQ ))  ~V{q), 

As  above,  our  nonholonomic  constraints  are  given  by  a  distribution  V  C  TQ.  We  also  let 
V°  C  T*Q  denote  the  annihilate  of  this  distribution.  Using  a  basis  ua  of  the  annihilate  T>°, 
we  can  write  the  constraints  as 

wa(q)  =  0, 


where  a  =  1, . . . ,  k. 

Recall  that  the  cotangent  bundle  T*Q  is  equipped  with  a  canonical  Poisson  bracket  and  is 
expressed  in  the  canonical  coordinates  (q,p)  as 


{F,G}{q,p) 


dFdG  dFdG 

dq1  dpi  dpi  d(p 


Here  J  is  the  canonical  Poisson  tensor 


fdFP  dFp\  T(  f 

\  dq  ’  dp  ) 


(  O  n  In  \ 
\  In  On  J 


A  constrained  t  phase  space  At  —  F L(T>)  C  T*Q  is  defined  so  that  the  constraints  on  the 
Hamiltonian  side  are  given  by  p  E  At.  In  local  coordinates, 

M  =  l(q,p)eT*Q\ojt^-  =  o\. 

Let  {Xa}  be  a  local  basis  for  the  constraint  distribution  V  and  let  {uja}  be  a  local  basis  for 
the  annihilator  V°.  Let  {c ua}  span  the  complementary  subspace  to  V  such  that  (uja1  uji)  =  5%, 
where  5^  is  the  usual  Kronecker  delta.  Here  a  —  1, . . . ,  k  and  a  =  1, . . . ,  n  —  k.  Define  a 
coordinate  transformation  (q,p)  i— >■  {q.  Pen  Pa)  by 

Pa  =  XlaPi ,  pa  =  UpjPi- 

In  the  new  (generally  not  canonical)  coordinates  (qiPonPa),  the  Poisson  tensor  becomes 

=)=  ( !'/•'/  !  {q\Pj} 

V I  v  •  ■/  }  !/'•/'  ! 


Use  (q,pa)  as  induced  local  coordinates  for  At.  It  is  easy  to  show  that 


dH  dHM 

AQiPcnPa)  =  -7p—(qiPa), 


dqJ 

dH 


dq] 

(q,Pa,Pa)  =  Hp(q,Pa), 


where  Hm  is  the  constrained  Hamiltonian  on  At  expressed  in  the  induced  coordinates.  We 
can  also  truncate  the  Poisson  tensor  J  by  leaving  out  its  last  k  columns  and  last  k  rows  and 
then  describe  the  constrained  dynamics  on  A4  expressed  in  the  induced  coordinates  (. q\pa )  as 
follows: 


^Adiq^  Pa) 


G  M. 


Here  Jm  is  the  (2 n  —  k)  x  (2 n  —  k)  truncated  matrix  of  J  restricted  to  Ai  and  is  expressed 
in  the  induced  coordinates. 


The  matrix  J ^  defines  a  bracket  {•  ,  on  the  constraint  submanifold  A4  as  follows: 

{Fm,Gm}m(Q,P„)  :=  -W,Pa)  (  S  )  > 

for  any  two  smooth  functions  F A4,  Gm  on  the  constraint  submanifold  A4.  Clearly,  this  bracket 
satisfies  the  first  two  defining  properties  of  a  Poisson  bracket,  namely,  skew  symmetry  and  the 
Leibniz  rule,  and  one  can  show  that  it  satisfies  the  Jacobi  identity  if  and  only  if  the  constraints 
are  holonomic.  Furthermore,  the  constrained  Hamiltonian  II m  is  an  integral  of  motion  for  the 
constrained  dynamics  on  A4  due  to  the  skew  symmetry  of  the  bracket. 


Following  e.g.  van  der  Sclmft  and  Maschke  [1994]  and  Koon  and  Marsden  [1997]  we  can  write 
the  nonholonomic  equations  of  motion  as  follows: 
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Jacobiator  of  the  Poisson  tensor  vanishes  precisely  when  the  curvature  of  the  nonholonomic 
constraint  distribution  is  zero  or  the  constraints  are  holonomic. 


•Euler-Poincare-Suslov  Equations 

Important  special  case  of  the  reduced  nonholonomic  equations. 

•Example:  Euler-Poincare-Suslov  Problem  on  50(3)  In  this  case  the  problem  can 
be  formulated  as  the  standard  Euler  equations 

Iuj  =  Iuj  X  UJ 

where  uj  =  (ctq,  cj2,  uj^)  are  the  system  angular  velocities  in  a  frame  where  the  inertia  matrix  is 
of  the  form  I  =  diag(/i,  /2, 13)  and  the  system  is  subject  to  the  constraint 

a  •  uj  =  0 


where  a  =  (flq,  a2,  as). 

The  nonholonomic  equations  of  motion  are  then  given  by 

Iuj  =  Iuj  x  uj  +  Xa 

subject  to  the  constraint.  Solve  for  A: 

I~la  •  {Iuj  x  uj) 

I~la  •  a 

If  a  is  an  eigenvector  of  the  moment  of  inertia  tensor  flow  is  measure  preserving. 


Invariant  Measures  of  the  Euler-Poincare-Suslov  Equations  An  important  spe¬ 
cial  case  of  the  reduced  nonholonomic  equations  is  the  case  when  there  is  no  shape  space  at  all. 
In  this  case  the  system  is  characterized  by  the  Lagrangian  L  —  abGaGb  and  the  left-invariant 
constraint 

(a,  tt)  =  aAttA  =  0.  (0.20) 

Here  a  =  a,AeA  £  0*  and  Q  =  (l  'c  i,  where  <■  _ \ ,  A  =  1, . . . ,  k,  is  a  basis  for  0  and  eA  is  its 
dual  basis.  Multiple  constraints  may  be  imposed  as  well.  The  two  classical  examples  of  such 
systems  are  the  Chaplygin  Sleigh  and  the  Suslov  problem.  These  problems  were  introduced 
by  Chaplygin  in  1895  and  Suslov  in  1902,  respectively. 

We  can  consider  the  problem  of  when  such  systems  exhibit  asymptotic  behavior.  Following 
Kozlov  [1988]  it  is  convenient  to  consider  the  unconstrained  case  first.  In  the  absence  of 
constraints  the  dynamics  is  governed  by  the  basic  Euler- Poincare  equations 

pB  =  Ccab\adpcPd  =  CcABpcCA  (0.21) 

where  ps  =  IabCb  are  the  components  of  the  momentum  p  £  gt  One  considers  the  question  of 
whether  the  (unconstrained)  equations  (0.21)  have  an  absolutely  continuous  integral  invariant 
fdkC  with  summable  density  Ai.  If  Ad  is  a  positive  function  of  class  C 1  one  calls  the  integral 
invariant  an  invariant  measure.  Kozlov  [1988]  shows 

Theorem  0.4  The  Euler- Poincare  equations  have  an  invariant  measure  if  and  only  if 
the  group  G  is  unimodular. 

A  group  is  said  to  be  unimodular  if  it  has  a  bilaterally  invariant  measure.  A  criterion 


for  unimodularity  is  CAC  =  0  (using  the  Einstein  summation  convention).  Now  we  know 
(Liouville’s  theorem)  that  the  flow  of  a  vector  differential  equation  x  =  f(pc)  is  phase  volume 
preserving  if  and  only  if  div  /  =  0.  In  this  case  the  divergence  of  the  right  hand  side  of  equation 
(0.21)  is  CACIADpB  —  0.  The  statement  of  the  theorem  now  follows  from  the  following  theorem 
of  Kozlov  [1998]:  A  flow  due  to  a  homogeneous  vector  field  in  W1  is  measure-preserving  if 
and  only  if  this  flow  preserves  the  standard  volume  in  Mn. 

Now,  turning  to  the  case  where  we  have  the  constraint  (0.20)  we  obtain  the  Euler- Poincare- 
Suslov  equations 

pB  =  CABlADpcPD  +  A  aB  =  CcABpcfilA  +  \aB  (0.22) 

together  with  the  constraint  (0.20).  Here  A  is  the  Lagrange  multiplier.  This  defines  a  system  on 
the  subspace  of  the  dual  Lie  algebra  defined  by  the  constraint.  Since  the  constraint  is  assumed 
to  be  nonholonomic,  this  subspace  is  not  a  subalgebra.  One  can  then  formulate  a  condition  for 
the  existence  of  an  invariant  measure  of  the  Euler-Poincare-Suslov  equations. 

Theorem  0.5  Equations  (0.22)  have  an  invariant  measure  if  and  only  if 

Kdudfliaa  +  T  =  / m ,  /iGl,  (0.23) 

where  K  —  l/(a,I_1a)  and  T  E  Q*  is  defined  by  (T,  f)  =  Tr(ad^). 

This  theorem  was  proved  by  Kozlov  [1988]  for  compact  algebras  and  for  arbitrary  algebras  by 
Jovanovic  [1998].  In  coordinates,  condition  (0.23)  becomes 

KCABfADacaD  +  CBC  =  paB. 


For  a  compact  algebra  (0.23)  becomes 

[I-1a,  a]  =  pa ,  p  G  M,  (0.24) 

where  we  identified  g*  with  g. 

The  proof  of  theorem  0.5  reduces  to  the  computation  of  the  divergence  of  the  vector  field  in 

(0.22). 

In  the  compact  case  only  constraint  vectors  a  which  commute  with  I ~la  allow  the  measure  to 
be  preserved.  This  means  that  a  and  I-1  a  must  he  in  the  same  maximal  commuting  subalgebra. 
In  particular,  if  a  is  an  eigenstate  of  the  inertia  tensor,  the  reduced  phase  volume  is  preserved. 
When  the  maximal  commuting  subalgebra  is  one-dimensional  this  is  a  necessary  condition. 
This  is  the  case  for  groups  such  as  50(3). 

We  thus  have  the  following  result  which  reflects  a  symmetry  requirement  on  the  constraints: 

Theorem  0.6  A  compact  Euler- Poincare- Suslov  system  is  measure  preserving  (i.e.  does 
not  exhibit  asymptotic  dynamics)  if  the  constraint  vectors  a  are  eigenvectors  of  the  inertia 
tensor,  or  if  the  constrained  system  is  Z2  symmetric  about  each  of  its  principal  axes.  If 
the  maximal  commuting  subalgebra  is  one- dimensional  this  condition  is  necessary. 


Invariant  Measures  of  Systems  with  Internal  Degrees  of  Freedom  In  this 
section  we  extend  the  result  of  Kozlov  [1988]  and  Jovanovic  [1998]  to  nonholonomic  systems 
with  nontrivial  shape  space.  One  can  think  of  these  systems  as  the  Euler- Poincare-Suslov 
systems  with  internal  degrees  of  freedom.  Recall  that  the  constraints  are  of  the  form  Qm+1  = 
•  •  •  =  Qk  =  0.  To  simplify  the  exposition,  we  consider  below  systems  with  a  single  constraint. 
The  results  are  valid  for  systems  with  multiple  constraints  as  well. 

Consider  a  nonholonomic  system  with  the  reduced  Lagrangian  Z(r,  r,  0)  and  a  constraint 
(a(r),  0)  =  0.  The  subspace  of  the  Lie  algebra  defined  by  the  constraint  at  the  configuration  q 
is  denoted  here  by  Qq.  The  orientation  of  this  subspace  in  q  depends  on  the  shape  configuration 
of  the  system,  r.  The  dimension  of  $q  however  stays  the  same.  As  discussed  in  section  ,  we 
choose  a  special  moving  frame  in  which  is  spanned  by  the  vectors  e\ (r), . . . ,  e^-i(r).  The 
equation  of  the  constraint  in  this  basis  becomes  Qk  =  0.  Recall  that  the  horizontal  part  of  the 
kinetic  energy  metric  is  gap(r). 

Theorem  0.7  The  system  associated  with  the  reduced  Lagrangian  Z(r,  r,  O)  and  the  con¬ 
straint  (a(r),C)  =  0  has  an  integral  invariant  with  a  C 1  density  A4(r)  if  and  only  if 

f  Ika  \ 

(i)  [CL  -  ctwj  - gaSViaS  =  0, 

(ii)  the  form  \Vlf3  —  ga6\^Bka^  dr13  is  exact. 


Systems  with  One-Dimensional  Shape  Space.  Assume  that  condition  (i)  of  theorem 
0.7  is  satisfied.  In  this  case  the  equation  for  the  density  of  the  invariant  measure  becomes 

d(\n  M)  =  d(\ng)  +  Vbbdr.  (0.25) 

The  solution  of  this  equation  is  globally  defined  if  the  shape  space  is  either  noncompact  (and 
thus  diffeomorphic  to  M),  or  compact  and  the  average  of  the  function  T>b  equals  zero. 

Systems  with  Conserved  Momentum.  If  the  nonholonomic  momentum  is  a  constant 
of  motion,  then  condition  (i)  of  theorem  0.7  is  trivially  satisfied.  Moreover,  condition  (ii)  now 
asks  that  the  form 

gai\uBka/ 

is  exact.  The  system  thus  preserves  the  measure  with  the  density 


(0.26) 


Examples 

The  Routh  Problem.  This  mechanical  system  consists  of  a  uniform  sphere  rolling  without 
slipping  on  the  inner  surface  of  a  vertically  oriented  surface  of  revolution.  He  described  the 
family  of  stationary  periodic  motions  and  obtained  a  necessary  condition  for  stability  of  these 
motions.  Routh  noticed  as  well  that  integration  of  the  equations  of  motion  may  be  reduced 
to  integration  of  a  system  of  two  linear  differential  equations  with  variable  coefficients  and 
considered  a  few  cases  when  the  equations  of  motion  can  be  solved  by  quadratures.  Modern 
references  that  treat  this  system  are  Hermans  [1995]  and  Zenkov  [1995]. 

This  problem  is  50(2)  x  SO (2)-invariant,  where  the  first  copy  of  50(2)  represents  rotations 
about  the  axis  of  the  surface  of  revolution  while  the  second  copy  of  50(2)  represents  rotations 
of  the  sphere  about  its  radius  through  the  contact  point  of  the  surface  and  the  sphere. 


Let  r  be  the  latitude  of  this  contact  point,  a  be  the  radius  of  the  sphere,  c(r)  +  a  be  the 
reciprocal  of  the  curvature  of  the  meridian  of  the  surface,  and  b(r)  be  the  distance  from  the 
axis  of  the  surface  to  the  sphere’s  center.  The  shape  metric  is  c2(r)r2/ 2  while  the  momentum 
equations  are 

c(r)  sin  r  2  f  c(r)  cos  r\ 

Pl  =  b(r)  PlT  ~  7P2r’  P2  =  V  “  W)  ) PlT' 

The  shape  space  is  one-dimensional,  the  symmetry  group  50(2)  x  50(2)  is  commutative, 
and  there  are  no  terms  proportional  to  r2  in  the  momentum  equations.  The  trace  term  in  (0.25) 
equals  c(r)  sinr/6(r),  and  thus  the  density  of  the  invariant  measure  for  the  Routh  problem  is 

M  =  <?{r)e  ~a7Tdr .  (0.27) 

The  group  action  in  this  problem  is  singular:  the  intersection  points  of  the  surface  of  revolu¬ 
tion  and  its  axis  have  nontrivial  isotropy  subgroups.  The  shape  coordinate  r  equals  ±7t/2  at 
these  points.  As  a  result, 

lim  A4(r)  =  lim  A4(r)  =  oo. 

— 7r/ 2  r—>  7t/2 


r 


The  Falling  Disk.  Consider  a  homogeneous  disk  rolling  without  sliding  on  a  horizontal 
plane.  This  mechanical  system  is  50(2)  x  5E?(2)-invariant;  the  group  50(2)  represents  the 
symmetry  of  the  disk  while  the  group  SE( 2)  represents  the  Euclidean  symmetry  of  the  overall 
system. 

Classical  references  for  the  rolling  disk  are  Vierkandt  [1892],  Korteweg  [1899],  and  Appel 
[1900].  In  particular,  Vierkandt  showed  that  on  the  reduced  space  V/SE( 2)-  the  constrained 
velocity  phase  space  modulo  the  action  of  the  Euclidean  group  SE( 2) — most  orbits  of  the 
system  are  periodic. 

The  shape  of  the  system  is  specified  by  a  single  coordinate — the  tilt  of  the  disk  denoted  here 
by  r.  The  momentum  equations  are 


Pi 

P2 


cos  r 


+ 


sin2  r 


mR 2  +  B  A  cos  r 


P2  r, 


Acosr 


sm  r 

Pi  +  -7 - P2 

Acosr 


r. 


Hence,  the  trace  terms  Vhh  in  (0.25)  vanish,  and  the  density  of  the  invariant  measure  equals  the 
component  of  the  shape  metric  g(r).  The  latter  equals  the  moment  of  inertia  of  the  disk  with 
respect  to  the  line  through  the  rim  of  the  disk  and  parallel  to  its  diameter.  Since  the  density 
of  the  measure  is  determined  up  to  a  constant  factor,  we  conclude  that  the  dynamics  preserves 
the  reduced  phase  space  volume. 


The  3D  Chaplygin  Sleigh  with  an  Oscillating  Mass.  The  three-dimensional  Chap¬ 
lygin  sleigh  is  a  free  rigid  body  subject  to  the  nonholonomic  constraint  v3  =  0,  where  v3  is  the 
third  component  of  the  (linear)  velocity  relative  to  the  body  frame.  The  Lagrangian  of  this 
system  is 

\m  [(V)2  +  (v2)2  +  (V)2]  +  t  [^(n1)2  +  I2(Q2)2  +  I3(tt3)2'  . 

In  this  formula  M  is  the  mass  of  the  body,  Ij  are  the  eigenvalues  of  its  inertia  tensor,  and 
(n1,  fi2,  fi3)  and  (v\v2,v3)  are  the  angular  and  linear  velocities  relative  to  the  body  frame. 
The  dynamics  of  this  system  is  discussed  in  Neimark  and  Fufaev  [1972]. 

We  couple  this  system  with  an  oscillator  moving  along  the  third  coordinate  axis  of  the  body 
frame.  The  mass  of  this  oscillator  is  m  and  the  displacement  from  the  origin  is  r.  To  keep 
the  notation  uniform  with  the  general  theory,  we  write  the  components  of  the  linear  velocity 
relative  to  the  body  frame  as  (hi4,  hi5,  hi6).  The  vector  (hi1,  hi2,  hi3,  hi4,  hi5,  hi6)  should  be  viewed 
as  an  element  of  the  Lie  algebra  se( 3).  The  Lagrangian  of  this  new  system  is 

L  =  \  [h{Q})2  +  j2(02)2  +  /3(03)2]  +  [(fi4)2  +  (U5)2  +  (fi6)2; 

Zj  Zj 

rn 

+  —  [(fl4  +  fi2r)2  +  (Q5  -  Q}rf  +  (Q6  +  rf]  -  U(r).  (0.28) 

2 

The  configuration  space  is  M  x  SE( 3),  and  the  system  is  invariant  under  the  left  action  of  SE( 3) 
on  the  second  factor.  We  have  not  specified  the  potential  energy  as  its  choice  does  not  affect 
the  existence  of  the  invariant  measure.  The  shape  space  is  just  the  first  factor  of  1  x  SE( 3) 
and  is  one  dimensional,  and  thus  the  above  theory  is  applicable.  To  show  the  existence  of  the 


invariant  measure,  we  note  the  following: 

1.  The  constrained  Lagrangian  does  not  contain  terms  that  simultaneously  depend  on  r  and 
pa.  The  constraint  is  Q6  =  0.  Therefore,  all  the  coefficients  of  the  nonholonomic  connection 
as  well  as  its  curvature  form  vanish.  This  implies  that  the  terms  T>aap  and  )Cn/p  vanish  too. 
The  differential  form  from  condition  (ii)  of  theorem  0.7  is  therefore  trivial. 

2.  The  moving  frame  is  r-independent.  Therefore  all  of  the  coefficients  7^Q,  are  trivial.  Con¬ 
dition  (i)  of  theorem  0.7  is  satisfied  since  the  group  SE( 3)  is  unimodular  and  e§  is  the 
eigenvector  of  the  inertia  tensor. 

3.  The  shape  metric  is  r-independent. 

The  system’s  dynamics  preserves  the  volume  in  the  reduced  phase  space. 

This  can  be  verified  by  a  straightforward  computation  of  the  divergence  of  the  vector  held 
that  defines  the  equations  of  motion: 

dUa 

r  = - 

dr  ’ 

Pi  =  —Q2p3  T  hl3p2  —  ra£25r, 
p2  =  —Q  Pi  T  ft  P3  T  mO  f, 

P3  =  -&P2  T  h>2pi  -  h>4_P5  T  h>5£>4, 

Q  9 

P4  =  Q  P5  —  mVt  r, 

Q  "I 

P5  =  —O.  P4  T  rriQ  r. 


Chaplygin  Sphere.  This  system  consists  of  a  sphere  rolling  without  slipping  on  a  hor¬ 
izontal  plane.  The  center  of  mass  of  this  sphere  is  at  the  geometric  center,  but  the  principal 
moments  of  inertia  are  distinct.  Chaplygin  [1903]  proved  integrability  of  this  problem.  Modern 
references  for  the  Chaplygin  sphere  include  Kozlov  [1985]  and  Schneider  [2002]. 

One  may  view  this  system  as  a  nonholonomic  version  of  the  Euler  top.  The  configuration 
space  is  diffeomorphic  to  50(3)  x  M2.  We  choose  the  Euler  angles  (0,  <p,  pj)  and  the  Cartesian 
coordinates  (x,  y )  as  the  configuration  parameters  of  the  Chaplygin  sphere.  The  Lagrangian 
and  constraints  written  in  these  coordinates  become 


L  =  —  (0  cos  pj  T  <p  sin  pj  sin  6)2  +  — 
2  2 


■6  sin  pj  T  (p  cos  pj  sin  0)A 


+  +  T os  02  +  yO2  +  y2) 


and 


x  —  6  sin  <p  +  p)  cos  <p  sin  6  =  0,  y  T  6  cos  <p  +  ip  sin  <p  sin  6  =  0, 

respectively. 

This  system  is  SE( 2)  invariant.  The  action  by  the  group  element  (ce,  a,  b)  on  the  configuration 
space  is  given  by 

(0,  ip,  (p ,  x ,  y )  i— >  (6,  ip,(p  T  a,  x  cos  a  —  y  sin  a  +  a,  x  sin  at  y  cos  a  +  b). 


The  shape  space  for  the  Chaplygin  sphere  is  diffeomorphic  to  the  two-dimensional  sphere.  The 
nonholonomic  momentum  map  has  just  one  component  and  is  moreover  preserved.  Straight¬ 
forward  computations  show  that  the  form  (0.26)  is  exact.  The  conditions  for  measure  existence 


are  therefore  satisfied.  The  density  of  the  invariant  measure  is  computed  in  overdetermined 
coordinates  in  Chaplygin  [1903]  (see  also  Kozlov  [1985]). 

The  invariant  manifolds  of  the  Chaplygin  sphere  are  two-dimensional  tori.  The  phase  flow 
on  these  tori  is  measure  preserving  and  thus  there  are  angle  variables  (x,  y )  on  each  torus  in 
which  the  flow  equations  become 

A  .  fi 

M(x,y)’  V  M(x,y)' 

See  Kolmogorov  [1953]  and  Kozlov  [1985]  for  details.  In  general,  these  equations  cannot  be 
rewritten  as 

x  =  A,  y  =  /a. 

The  flow  however  becomes  quasi-periodic  after  a  time  substitution  dt  =  AA(x,  y)dr  (see  Kozlov 
[1985]  for  details).  This  example  thus  shows  that  the  flow  on  the  nonholonomic  invariant  tori 
can  be  more  complicated  than  a  Hamiltonian  flow. 

It  follows  that  adding  a  symmetry  preserving  potential  to  the  Lagrangian  of  the  Chaplygin 
sphere  leaves  the  new  system  measure  preserving  with  the  same  measure  density.  This  was 
pointed  out  by  Kozlov  for  a  specific  potential  (see  Kozlov  [1985]  for  details). 


A  Spherical  Robot:  Controlled  Chaplygin’s  Sphere  and  Chaplygin’s  Top 

Robot  viewed  as  a  controlled  Chaplygin’s  sphere  or  Chaplygin’s  top. 


Figure  0.7:  Schematic  configuration  of  a  spherical  robot  in  a  uniform  gravitational  field  with  a  slider  and  a  rotor 

Assume  the  spherical  base  body  can  roll  without  sliding  on  a  horizontal  plane  in  a  uniform 
gravitational  held;  the  radius  of  the  ball  is  a.  See  Figure  0.7  for  its  schematic  configuration. 
Choose  a  base  body  coordinate  frame  with  the  origin  at  the  center  of  the  ball.  Let  x  E  M3 
denote  the  position  of  the  center  of  the  ball  in  the  inertial  frame,  and  let  R  E  SO (3)  represent 
the  base  body  attitude  that  maps  from  the  base  body  coordinate  frame  onto  the  inertial  frame. 
The  spherical  robot  is  controlled  by  internal  actuators.  Relative  motion  of  the  internal  actuators 
with  respect  to  the  base  body  is  described  by  generalized  shape  coordinates  r  E  QS)  where  Qs 
is  referred  to  as  the  shape  space.  Hence,  the  configuration  space  manifold  is  SO(3)xM3  x  Qs. 
Note  that  we  are  only  interested  in  (aq,^),  the  horizontal  position  of  the  center  of  the  ball. 
The  general  theory  of  developing  equations  of  motion  in  the  form  of  Euler- Poincare  is  given  in 
the  next  subsection. 


General  System 

Let  G  be  a  general  linear  matrix  Lie  group  and  Lb  be  a  vector  space;  let  S  denote  the  semidi- 
rect  product  G(S)bb.  The  shape  space  Qs  is  an  n-dimensional  Abelian  Lie  group  with  local 
coordinates  r  =  (ri, • •  •  ,  rn).  The  conhguration  manifold  is  Q  =  S  x  Qs.  We  assume  that 
G  x  Qs  forms  a  trivial  principal  fiber  bundle  and  that  for  every  q  =  (g.  r)  E  Q  where  g  E  G 
and  r  E  QSl  the  left  action  of  G  on  G  x  Qs  is  a  smooth  map  <L  :  G  x  (G  x  Qs)  ^(Gx  Qs) 
given  by  T/gy  =  (. hg ,  r)  for  an  arbitrary  h  E  G,  where  T  is  assumed  to  be  free  and  proper.  We 
assume  that  controls  act  on  the  shape  space  Qs  so  that  the  shape  is  fully  actuated. 


Let  the  reduced  Lagrangian  be 


l(Z,Y,r,r,r)  =  T(Z,Y,r,r)-V(r,r), 

where  f  =  g~lg  E  0,  T  =  ^_1ao  E  Lb  for  a  fixed  ao  E  Lb.  Let  y(£)  be  a  (smooth)  curve  in  Lb, 
Y  =  E  Lb.  Assume  the  reduced  kinetic  energy  T(£,  Y,  r,  r)  can  be  written  as 


TK,Y,r,r)  =  -^,Y,rjM(r)  I  Y 


where  M(r )  denotes  a  reduced  inertia  tensor  on  only.  Assume  the  reduced  constraint  can 
be  written  as  Y  —  f£(r),  where  C  :  Lb  — >  Lb  is  a  smooth  vector-valued  function.  Thus  the 
constrained  reduced  Lagrangian  becomes 


/cK,r,r,r)  =  TcK,r,r,r)-Y(r,r), 


where 


Here  M(T,  r)  defines  another  reduced  inertia  tensor  dependent  on  Qs  and  the  dynamic  param¬ 
eter  T  G  W,  which  is  referred  to  as  “advected  parameter”  in  the  literature.  Since  T  =  g~1ao 
evolves  on  G,  one  may  view  M(T,  r)  as  a  function  on  G  x  Qs. 


Equations  of  Motion  of  the  Spherical  Robot 

In  this  section,  we  apply  the  general  results  given  in  the  previous  section  to  the  spherical 
robot.  For  this  example,  the  Lie  group  G  =SO(3)  and  the  vector  space  W  =  M3. 


Let  mo  be  the  mass  of  the  base  body,  Jo  be  the  inertia  tensor  of  the  base  body  defined  with 
respect  to  the  base  body  coordinate  frame,  rrq,  i  =  1,  •  •  •  ,  AT,  be  the  mass  of  the  z-th  auxiliary 
body,  and  Jt{r)  be  the  inertia  tensor  of  the  z-th  auxiliary  body  defined  with  respect  to  the  base 
body  coordinate  frame.  Moreover,  let  po  denote  the  relative  position  vector  of  the  center  of 
mass  of  the  ball  and  let  pi(r)  denote  the  relative  position  vector  of  the  center  of  mass  of  the 
z-th  auxiliary  body.  We  denote  the  angular  velocity  of  the  z-th  body  relative  to  the  base  body 
coordinate  frame  by  Ci(r)r.  That  is,  suppose  the  orientation  of  the  z-th  body  in  the  base  body 
frame  is  given  by  Ri(r)  E  50(3),  where  Ri(r)  maps  from  a  coordinate  frame  for  the  z-th  body 
into  the  base  body  coordinate  frame.  Let  ay(r,  f)  =  R~[l{r)j-fRi(r)  denote  the  angular  velocity 
of  the  z-th  body  relative  to  its  own  coordinate  frame.  Then  Ci(r)r  =  Ri(r)cJi(r ,  r). 

We  obtain  the  reduced  equations  of  motion  for  50(3)  x  Qs 


d  /dlc\  dlc 
dt\d(j)  Duj 

_  dlc 

dt\drJ  dr 


— ^ u  p  _  (pc(r)  x  T)  •  uj  + 
+mTag r  x  pc(r), 


(0.30) 


where  the  shape  dynamics  are  completely  controlled. 


